The Glass Transition in Driven Granular Fluids: A Mode-Coupling Approach 
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We consider the stationary state of a fluid comprised of inelastic hard spheres or disks under 
the influence of a random, momentum-conserving external force. Starting from the microscopic 
description of the dynamics, we derive a nonlinear equation of motion for the coherent scattering 
function in two and three space dimensions. A glass transition is observed for all coefficients of 
restitution, e, at a critical packing fraction, tp^[e), below random close packing. The divergence of 
timescales at the glass-transition implies a dependence on compression rate upon further increase 
of the density - similar to the cooling rate dependence of a thermal glass. The critical dynamics 
for coherent motion as well as tagged particle dynamics is analyzed and shown to be non-universal 
with exponents depending on space dimension and degree of dissipation. 

PACS numbers: 64,70Q-,81.05Rm 



I. INTRODUCTION 

A wide range of fluids can be quenched into a disor- 
dered, solid state. This includes metallic melts [1], col- 
loidal suspensions [2], foams [3] and recently, evidence 
was given that granular fluids may also undergo a glass 
transition [4-10]. Among all these different systems, col- 
loidal suspensions are probably best understood. Exper- 
iments by van Megen et al. [2, 11] showed that besides 
the fluid and the ordered crystalline phase, colloidal sus- 
pensions in thermal equilibrium can also form colloidal 
glasses: A dynamically arrested state of the system which 
is characterized by diverging relaxation times [12]. While 
a complete theoretical understanding of the glass transi- 
tion in fragile glass formers is still missing [13], mode 
coupling theories can quite successfully describe some of 
the phenomena on a quantitative level [14]. 

One interesting question has been raised more recently: 
Docs the glass transition survive, if the system is driven 
by external forcing into a nonequilibrium state? Or 
more generally, can one observe a glass transition also 
in a nonequilibrium system? A well studied example are 
sheared colloidal suspensions for which it was shown that 
the equilibrium glass transition disappears [15, 16]. An- 
other example is nonlinear microrheology [17-19], where 
a strong pulling force is applied to a single particle, forc- 
ing it out of its cage, thereby possibly melting the glass. 

Another system far from equilibrium are athermal 
packings of particles [20, 21], undergoing a jamming tran- 
sition at a critical packing fraction. Many of the proper- 
ties close to the jamming point resemble those of fluids 
at the glass transitions. This observation is at the heart 
of the jamming diagram, where the glass transition in 
thermal systems and the jamming transition are part of 
a larger parameter space [22, 23]. 

Since granular particles are too large to be thermally 
activated, one necessarily needs a driving force to keep 
the grains in motion for extended periods of time. While 



in nature, gravity is probably the most important driv- 
ing force [24], experimentalists have devised quite a few 
methods of fluidisation. The list includes shaking [25], 
electrostatic [26, 27] or magnetic [27, 28] excitation and 
fluidization by air [29, 30] or water [31]. 

We have recently investigated the possibility of a glass 
transition in driven granular fluids. In two publications 
[9, 10], henceforth referred to as I and II, we have demon- 
strated that mode coupling theory (MCT) can be general- 
ized to the far from equilibrium stationary state of a gran- 
ular fluid. In particular we found a granular glass tran- 
sition for all degrees of dissipation, accompanied by the 
common signatures of a dense fluid close to the glass tran- 
sition. Here, a careful derivation of the granular mode 
coupling equations is presented and the consequences 
of MCT are worked out in detail. We furthermore ex- 
tend our previous analysis to two-dimensional systems, 
which are realized in many experiments on granular mat- 
ter [5, 7, 8, 25, 30]. The resulting glass transition diagram 
is shown in Fig. 1 in the plane spanned by packing frac- 
tion, if, and coefficient of restitution, e. 

The dissipative interactions of the granular particles 
imply two primary consequences. First, while the dynam- 
ics of particles in thermal equilibrium is microscopically 
time reversal invariant, the symmetry under time rever- 
sal is broken for granular dynamics. Second, there is no 
natural equilibrium reference state like for the sheared 
colloids [15, 16], were the fluid can be thought of as be- 
ing driven out of equilibrium by the optional external 
driving force. In the granular system, the driving force is 
required to maintain a stationary state with more than 
transient dynamics. 

The paper is organized as follows. In Sec. II we 
deflne the model of a driven granular fluid and intro- 
duce the microscopic dynamics in Sec. III. In Sec. IV 
we derive the MCT equations for the coherent scatter- 
ing function, <j){q,t), discuss the asymptotic correlations 
fq :~ (j){q, t oo), used as an order parameter to locate 
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FIG. 1. Critical packing fraction, tpc, separating the fluid from 
the glassy state of driven granular particles as a function of 
the coefficient of restitution e for space dimension D = 2 (top) 
and D — 3 (bottom). 

the glass transition, and analyze the dynamics close to 
the glass transition. The MCT equations for the incoher- 
ent scattering function, (j)s{q,t), of a tagged particle are 
derived in Sec. V. In Sec. VI we discuss our results in a 
broader context and conclude with a number of perspec- 
tives for future work in Sec. VII. 



II. MODEL 

A. Inelastic Hard Spheres 

The granular fluid is modeled as a monodispersc sys- 
tem of N smooth inelastic hard spheres (in dimension 
D ~ 3) or disks (in D = 2) of radius a and mass m ~ 1 
in a volume V = . We consider the thermodynamic 
limit N,V ^ oo such that the density n = N/V remains 
finite. Dissipation is introduced through a constant coef- 
ficient of normal restitution eg [0; 1] that augments the 
law of rcfiection [32], 

1^12 • ■y'l2 = -£'^12 • ■"12, (1) 

where Vi2 = V1—V2 is the relative velocity and ri2 is the 
unit vector pointing from the center of particle 2 to par- 
ticle 1. The prime indicates post-collisional quantities. 

B. Stochastic Driving Force 

The driving force is implemented as an external ran- 
dom force, 

v'.it) = v,it) + ^Ut). (2) 



where Pd is the driving power. The , a = 1, . . . , D are 
Gaussian random variables with zero mean and variance, 

(i'))^ = - 5^i.v>,,]5^^5{t - t'), (3) 

where 7r(i) = argminj.{|ri — r^l > £d} yields the index of 
the particle that is closest to particle i but at least a given 
distance, £d, away. Thereby, the external force does not 
destroy momentum conservation on length scales i ^ £d- 
We choose in on the order of a mean particle separation. 

C. The Granular Fluid 

For undriven granular fiuids, it is known that the homo- 
geneous cooling state is unstable to shear- and eventually 
to density fiuctuations [33, 34]. In fact, the particles form 
extremely dense clusters. No such clustering instability 
is predicted [33], and indeed observed, for the randomly 
driven fluid. Consequently, we assume that the fluid is 
macroscopically homogeneous and isotropic. This implies 
that all spatial two-point correlation functions C{r,r') 
are functions of the distance |i — r'\ only. In the station- 
ary state, the system is also time translation invariant, 
implying that time dependent correlation functions are 
only functions of time differences. 

Macroscopically, the fluid is fully characterized by the 
packing fraction, ip (where ip = Anna'^ /S in Z? = 3 and 
ip = irna^ in D = 2), the coefficient of restitution, e, and 
the driving power, Pp. A more conventional description 
is given in terms of the packing fraction and the granular 
temperature T = T{(p,e, Pd) = jjj^J^i'^i- The latter 
is given by the balance between the driving power, Pq, 
and the energy loss through the inelastic collisions, P^ :~ 
U,{^,e)u;,T, where U,{^,e) ~ U.^sie) = (1 - e^)/D in 
a mean field approximation [32] . 

The collision frequency ojc oc VT is the only time scale 
of the system. Thus, changing the granular temperature 
only changes the time scale of the system. To keep the 
discussion more transparent, we refrain from using the 
freedom to set T = 1 but keep in mind that the qualita- 
tive behavior of the system is independent of the temper- 
ature T > 0. 



III. MICROSCOPIC DESCRIPTION 

A. Observables 

The two relevant observables discussed in the follow- 
ing are the density field, p(r,i), and the current density, 
j{r,t), with the following microscopic definitions: 

i 

3{r,t) = ^Y.''^mr~r,{t)). (4b) 
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We will use the spatial Fourier transforms Pk{t) = 
FT[p]{t) and the longitudinal current j^{t) = k ■ FT[j](t). 
The corresponding tagged particle quantities are given as 

p%r,t)=6{r-rs{t)), (5a) 
f{r,t)=vj{r-r,{t)). (5b) 

B. Dynamics 

The (forward in time) pseudo Liouville operator de- 
scribes the time evolution of a microscopic observable A, 
i.e., iL^A = dA/dt, according to the dynamics specified 
above [35]. It is given as the sum of three parts, 

£+=£o+^T+ +£+, (6) 
j<fe 

which are in turn: (i) The free streaming operator iLq — 
"Y^j Vj ■ V j. (ii) The collision operator, 

= -(^jfc ■ VjkjQi-fjk ■ Vjk)5(rjk - 2a)(6+ - 1), (7) 

where Q{x) denotes the Heavisidc step function and the 
operator implements the inelastic collision rule [36] 
and (iii) the driving operator [37], 




With the binary collision expansion [38] , formal power 
series of the pseudo Liouville operator can be de- 
fined. In particular, this allows to write the propagator 
U(t) = exp(zt£_|_) in terms of an exponential operator 
exp(it£+) := X]n(*^'^-i-)"/"-'- The Laplace transformed 
propagator U(s) = LT[U](s) = (s — £ + )^^ is then also 
defined as a power series [39] . 

The starting point for the derivation of equations of 
motion is an operator identity that is most concisely ex- 
pressed in the Laplace domain, 

TU(s)T= [s-Q- M(s)]-\ (9) 

where Q = VL+y, 

U{t) = T£+Qexp(itQ£+Q)Q£+T, (10) 

and y = y'^ ,0. — 1 — y are projection operators [40]. 

C. Phase Space Averages 

Our approach starts from the microscopic description 
of the particle dynamics in terms of the location in 
phase space, F := (ri, t^i, . . . rjv, i)jv) and the trajec- 
tory of the external driving force = (^i(t), . . . , ^jv(^))- 



Macroscopic observables A{r,t\'E.t) := (a(r, t|F, S())p = 
/ o?Fp(F, t)a(r, <|F, Sj) are then introduced as expecta- 
tion values with respect to the phase space distribution 
function p(F, t) of the microscopic variables. Here, we 
restrict ourselves to the stationary state, where the dis- 
tribution function is time independent, p(F,t) = g(r). 

In contrast to fluids in equilibrium, no analytical ex- 
pression for the stationary phase space distribution of 
driven granular fluids is known so far. Therefore we have 
to make a few assumptions to evaluate the expectation 
values. First of all we assume that positions and veloc- 
ities are uncorrelated, q(T) = Qr{{ri})Qv[{vi})- More- 
over, we assume that the velocity distribution factor- 
izes into a product of one particle distribution functions, 
Qv{{vi]) — Y\iQi{vi)- All we need to know about Qi{v) 
is that it has a vanishing first moment, / d^v vqi{v) ~ 
and a finite second moment, J d^vv'^gi{v) = DT < oo. 
The spatial distribution function, gi,.({ri}), enters the the- 
ory via static correlation function, as will be discusssed 
below. 

A common expectation in disordered systems is, that 
the macroscopic expectation values should be self averag- 
ing, i.e., independent of any specific disorder realization 
[41]. In our model, this applies to the stochastic driv- 
ing force St, i.e., we define macroscopic observables as 
averages over all realizations of the driving, A{r,t) = 

{a{r,t)) := (^A{r , t\Et)) ^ = J dEtP{Et)A{r,t\Et). Here, 

P{Et) is the distribution of the random forces. Aver- 
ages over pairs of observables define a scalar product, 
{A\B) := {A*B) where A* denotes the complex conju- 
gate of A. 

Given the definition of a scaler product, we can for- 
mally introduce the adjoint Liouville operator, /C]|_, via 

the relation ^£^a]6^ = {a\Z+b). For elastic hard spheres 
in thermal equilibrium, it can be shown that detailed 
balance implies £^(e = 1) = £_(e = 1), where is 
the Liouville operator describing time reversed dynamics. 
This relation does not hold for inelastic collisions which 
arc discussed here. In the present context, an explicit 
expression for £^ is not needed and hence will be given 
elsewhere. 

As both the hard sphere interactions and the 
driving force (by construction) conserve momentum, 
this will also be refiected in the matrix elements 
(ai(fci)---a„(fc„)16i(pi)---6,„(p,„)) Sj^k^.Y^P. the 
form of a selection rule. 

The central quantities in the following will be the co- 
herent scattering function 

0(g,i) ■.= N{pq{T)\pg{T + t))/S,, (11a) 

where Sq :~ N {pq\pq) is the static structure factor, and 
the incoherent scattering function 

M<l,t)--^{Pg{r)\p'qir + t)). (lib) 

In general, all macroscopic quantities will be functions of 
the coefficient of restitution e. To reduce clutter, we will 
suppress this dependence. 
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D. Static Structure Factors 

Below, we will treat the static structure functions as a 
known input. Hence, we need Sg = Sq{ip, e) for a range of 
densities, if around where we expect the critical density 
to be. Lacking good quality data for Sq{Lp ^ (pc{e),e), 
let alone reliable theoretical predictions for this quantity, 
we use preliminary results that Sq{ip,e) only weakly de- 
pends on the coefficient of restitution e and approximate 
Sq{(p,e) « Sq{ip,s = 1) by their elastic counterparts: In 
D = 3 we use the Percus-Yevick (PY) equation [42] for 
elastic hard spheres in thermal equilibrium [43], except 
for the pair correlation function at contact which is bet- 
ter approximated by the Carnahan-Starling expression 
[44]. In = 2 we use the Baus-Colot (BC) equation [45] 
throughout. 



where = fip^Tij-p, or, equivalently, in the time domain 
as the solution of the equation of motion. 



4>iq,t) + iyq4>{q,t)+nl(f>{q,t) 

+ 1 dTm{q,t-T)<p{q,T) 



+ / dTL{q,t~T)^{q,T) = 0, 



(16) 



where m(g, t) = M(g, t)/^q and the initial conditions are 

(}){q, t = 0) = and (j){q, t = 0) = 1. To proceed, we need 
to find approximate expressions for the memory kernels. 
Before we come to the mode coupling approximation, we 
discuss the simpler assumption that M{q,t) = L{q,t) = 
0. 



IV. THE GRANULAR GLASS TRANSITION 



A. Equations of Motion 



Let us introduce the following microscopic state vector 



iq = VN{Pq/ \/S^, jq /VT). Then the coherent scatter- 
ing function (f>{q, t) is given as one element of the matrix 
of correlators $((?, t) = (aq|U(t)aq). 

With the help of Eq. (9) and the projectors 

J'c = NY,\pq) {p^\ /Sq +Nj2\jq) {jq\ /T, (12) 



Qc = 1 — CPc, one finds 



s ~n„{q)[l+^L{q,s)] 
-^jp{q) s - Uq- M{q,s) 



(13) 



M{q,t)^N{FimFq)/T, 
L{q,t)^ljlmFq) lq{j^\j^), 



while Vlpp (X {pq\L+pq) = q{pq\jq) = duc to par- 
ity. The other entries of the frequency matrix flab = 
{aq\Q.aq) are nonzero [Note that i^q = ^ljj{q)]. 
The memory kernels are formally given as 

(14a) 
(14b) 

where U{t) = exp{itQcL+Qc) is a modified propagator. 
There are two fluctuating forces, Fq = Qc£+jq, and 

~ Qc£j^^j^, and at this point we can not rule out 
that there is a nonzero fluctuating current, jt = Q^/C^pg, 
while Jq = Qc£+Pq = qQcjq = 0. In the elastic limit 
F^ = Fq and = Jq = holds and therefore L{q,t) 
vanishes. 

In the Laplace domain, the coherent scattering func- 
tion is thus given as 



4> ^{q^s) = s 



L{q,s)] 



M{q,s) 



(15) 



B. Sound Waves 

The linear equation of motion 

4>{q,t)+iyq<P{q,t)+nlcj){q,t)=0 (17) 

describes damped sound waves (j){q, t) = 
e^^o^/'^cosiCqqt). 

Sound damping duc to collisions, as described by the 
second term in Eq. (17), 

(18) 

can be evaluated in the Enskog approximation. The cal- 
culation for two dimensions is shown in appendix B 1 , 
yielding 

1 + e 

l^q = - 



-ujE[l + 2Ji;{2aq)] in £> = 2, (19a) 



where Jo{x) is the zeroth order Besscl function [46] and 
the double prime denotes the second derivative with re- 
spect to the argument. The result in three dimensions is 
known [47], 

uq = ^-^lue[1 + Sj'^{2aq)] in D = 3, (19b) 

where jo{x) — sin(x)/x is the zeroth order spherical 
Bessel function [48]. The Enskog collision frequency 
loe — 2^ D(tpxl2a)\jT j-K is given in terms of the con- 
tact value, X, of the pair correlation function [49]. 

The speed of sound is given in the long wavelength 
limit by C| = ^qjq^ = VlpjVljp/ q^ . One finds. 



n,p=N{j^\L+Pq)/^q 

qSu{,q)\jT/Sq 

with the longitudinal current correlator Su(q) 
N{jq\jq) IT. The calculation (cf. appendix B 2) of 



(20a) 



■ {Pq\^+Jq) 



q^/T/Sq ( ^ + 



(20b) 
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FIG. 2. Speed of sound, Cq, according to Eq. (21) as a func- 
tion of wave number, q, in 3D for packing fraction tp — 0.516 
and coefficient of restitution e — 1.0 (squares), 0.5 (filled cir- 
cles), and 0.0 (diamonds). 



uses the approximate granular Yvon-Born-Green (YBG) 
relation (of. appendix A). Combining these results, we 
find that the long wavelength speed of sound. 



c: = T 



Su{q ) fl + 

Sn 



2 



(21) 



is reduced for the dissipative driven fluid compared to a 
fluid of elastic hard spheres in thermal equilibrium. The 
sound damping, Vqjlq^ ^ on the other hand decreases with 
increasing dissipation. 

Alternatively, the speed of sound can be given as C| = 
S u{q) I nn'^^ ^ where the effective compressibility is 
defined in a form, 



1 



in D = 2. presumably implying logarithmic corrections 
to the sound damping in D = 2. 

In the following, we set Su{q) = 1. In Fig. 2, the 
resulting sound dispersion relations are shown for packing 
fraction = 0.516 and coefficient of restitution varying 
from £ = 1 to 0. The speed of sound decreases with 
increasing dissipation in agreement with hydrodynamie 
predictions [52]. 



C. The Mode Coupling Approximation 

In the spirit of the equilibrium mode coupling theories 
[53-56], we introduce a second projection operator 

^2 = ^'^|PfcPp)(PfcPp|/^fc^p (23) 

and approximate the modified propagator as 
U(t) « 72^)72 

where in the second step, a factorization approxima- 
tion, N'^ (^PkPp\y)it)pkPpj /SkSp « (j){k,t)4){p,t), was 

employed. Eq. (24) is known as the mode coupling ap- 
proximation (MCA). 

1. The MCA of L 



We find 



N 



Hq, i)-^J2 ^'jfcpW<jfcp0(A:, mip, i), (25) 



-Tnh 



(22) 



where 



reminiscent of a random phase approximation [50]. Here, 
= Sq/nT is the compressibility of a fictitious elastic 
hard sphere system with a structure factor Sq and Cq is 
the direct correlation function [49]. 

We close this section with a few remarks. First, in a 
thermal fiuid the expression for the sound velocity simpli- 
fies, because Su{q) = 1 due to molecular chaos [51]. In a 
granular fluid, Su{q) is actually found to be wave number 
dependent. Preliminary results indicate Su{q — >■ 0) < 1. 
Second, from ^ Qjp for e < 1, it can be seen ex- 
plicitly that the Liouville operator is not self adjoint. In 
terms of physical processes, this reflects that the transi- 
tion rate for the conversion of density fluctuations into 
current fluctuations is not equal to the rate of the reverse 
process. Detailed balance, or more general, time rever- 
sal invariance it broken already for the linear equation 
of motion (cf. Sec. VA below). Finally, it is known [49] 
that Navier-Stokes-order hydrodynamics does not exist 



Uqkp = N {pq\L+QcPkPp) /Sk^O 



(26) 



due to parity and Wqkp is defined below. Therefore 
L{q,t) = within the mode coupling approximation. 



2. The MCA of M 



The mode coupling approximation for M{q,t) yields 



N ^ 

M[q, t)^-Y, ^qkpy^qkpHk, t)(P{p, t), 



T 



k.p 



where 



Wqfcp = N {pkPp\QcJl^+]^) /Sp. 



(27) 



(28a) 
(28b) 
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The left vertex is known from the hterature [57] 
V„kp - -^[{q-k)Sp + {q-p)Sk-qS'-'Hk,p)/S,]Sg.k+p. 

Customarily, the convolution approximation [58], 
S'-^^k,p) « SkSpSq, is applied to yield 

T 

Vqfcp = J^^pHq ■ k)nck + (g • p)ncp]5qM+p- (29a) 

By a nontrivial calculation (see appendix B 3), we were 
able to show that the right vertex is approximately given 
as 

1 + e T 

Wqfcp w -^ — Sk[{q-k)nck + {q-p)ncp\5q^k+p, (29b) 
I 



different fromVqfcp. 



The physical interpretation of these results for the ver- 
tices is, that (i) the rate of annihilation of pairs of density 
fluctuations pk , Pp is determined by the static structure 
of the fluid, both in an equilibrium fluid and in the driven 
granular fluid; (ii) The rate of creation of such density 
fluctuations is suppressed by a factor (1 + e)/2, though, 
compared to the rate of creation or to the equivalent rate 
in an equilibrium fluid. 



The reduced memory kernel in the mode coupling ap- 
proximation (first reported in I) then reads 



m[<j>\{q,t) = A,{e)^ / d''kSkS\q_k\{[q-k]ck + [q-{q-k)]cpf4>{k,t)4>{\q-klt), (30) 



where 

A-\e) = l + \=^^S,. (31) 

Fig. 3 demonstrates the prefactor Aq{e) that distin- 
guishes the granular memory functions from the well- 
known elastic results where Aq{e = 1.0) = 1: For e < 1, 
the prefactor exhibits deviations from unity with oscilla- 
tions given by the static structure factor. As ^q(e) is 
minimal for the first peak of the structure factor, i.e. the 
length scale of the cage, one concludes that increasing 
dissipation (decreasing coefficient of restitution) weakens 
the cage effect. Compared to the elastic case, the force 
acting by the cage onto the particles inside the cage is 
smaller as the particles' reflections from each other are 
reduced by the influence of dissipation. While additional 
changes are expected by the e-dependence of the struc- 
ture factors, the major difference is encoded in the pref- 
actor Aq{e). The fact that Aq{£) > for all values of 
the coefficient of restitution e ensures that the memory 
kernel remains positive. 




FIG. 3. Prefactor of the memory kernel, Aq{e), as a func- 
tion of wave number, q, for two values of the coefficient of 
restitution, e = 0.5 (filled circles) and 0.0 (diamonds), in 3D. 



D. The approximate equation of motion and the 
phase diagram 

With the mode coupling approximations in place, the 
equation of motion 

4>iq,t) + iyq<Piq,t)+nlq^{q,t) 

/•* . (32) 

+ / drm[0](g,t-r)0(g,T) =0 
Jo 

turns into a closed equation for the coherent scattering 
function once the static structure factor, Sq, is known. 



This equation of motion has the same formal structure 
as the one for the elastic hard sphere fiuid in thermal 
equilibrium. The viscous term, Eq. (18), decreases with 
decreasing coefficient of restitution s; the speed of sound, 
Eq. (21), acquires a nontrivial dependence on the coeffi- 
cient of restitution as does the memory kernel m [(/)]. 

Structural arrest of the grains in a glassy state gives 
rise to time persistent density correlations. Hence, we 
introduce the order parameter for the glass transition, 
fq :~ limt_>.oo 0((7, t). It can readily be shown that the 
above equation of motion yields the following equation 
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for the asymptotic function, /g, 

T^-^ifM- (33) 

With the memory kernel being independent of tempera- 
ture, the order parameter fq is also independent of tem- 
perature as expected. It can easily be checked, that 
/q = is always a solution of the above equation. Study- 
ing this equation from a dynamical systems point of view, 
one finds that at a critical density ipc, the vanishing solu- 
tion becomes unstable and a new, stable solution /g > 
appears discontinuously, signaling the glass transition 
[14]. The order parameter at the critical density ipc will 
be denoted as f^. 

Using the structure factors as discussed in Sec. HID, 
we find the phase diagrams in Fig. 1. For D — 3 this 
result was first reported in I, whereas for D = 2 this is 
a new result (for technical parameters, cf. appendix C). 
The order parameter jumps discontinuously at the crit- 
ical density ipc as expected [14] from the type of singu- 
larity in Eq. (33). The evolution of the transition with 
e is remarkably similar for 2D and 3D, with transition 
densities increasing from the elastic case to e = 0.0 by 
around 10% in a roughly linear fashion. 




logt 



FIG. 4. Coherent scattering functions, (j>{Q, t), as a function of 
time, t, in 3D for wave number q — 4.2 /2a. At the transition 
point at packing fraction 'p'^{e — 0.5) = 0.548 with a critical 
glass-form factor of — 0.400, and at higher (1.1<^'^, 1.01^5'^, 
l.OOliys'^) and lower (0.9(y3^ O.QQvp^ O.QQQip") packing fractions. 
The exponent parameter is A = 0.710 which yields the critical 
exponent a = 0.323, and von Schweidler exponent b = 0.624. 
The critical amplitude for 2qa — 4.2 is hq = 0.583. The 
critical law labeled i"" is shown dashed for a time scale to = 
0.0260. 



E. Coherent Dynamics Close to the Glass 
Transition 

The full dynamics of density fluctuations is obtained 
by solving the MCT equations by iteration (for details, 
see Appendix C). In Figs. 4 and 5 we show the coherent 
scattering function in D = 3 for several densities, below 
and above the critical point for e = 0.5 and two wave 
numbers 2qa = 4.2 and 2qa = 7. As the critical point 
is approached from the fluid side, one observes the devel- 
opment of a plateau in the coherent scattering function. 
Increasing the density above the critical value (pc leads 
to an increase in the order parameter, fq [59]. 

The MCT equations of motion are known to admit 
scaling solutions at densities close to the critical density 
(pc{£ = 1) [14]. As the granular mode coupling equations 
are formally identical to those for an equilibrium hard 
sphere fluid, the scaling analysis is also closely related. 
To keep the presentation self contained, wc will summa- 
rize the main results of this analysis. The analogous anal- 
ysis for the incoherent scattering function discussed below 
in Sec. V was presented in II. 

At small distances a{e) := [fd^) — 'fi]/'Pci£) to the 
critical point, one finds 

Hq,t;a)=rq+hqG^it), (34) 

where hq = hq{e) is the critical amplitude and Gaif) is a 
scaling function, independent of the wave number. The 
scaling function Ga{t) can be characterized by a hierar- 
chy of time scales. The shortest time scale is naturally 
provided by the mean time between collisions, = u!~^. 
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FIG. 5. Coherent scattering functions, (j}{q,t), as a function 
of time, t, in 3D (cf. Fig. 4) for a different wave number 
q = 7/2a. The critical amplitude for 2qa = 7 is hq — 0.367. 
The von Schweidler law is shown dashed for a fitted time scale 
T = 1.21 X 10^ 

The time scale for the transition through the plateau, 
where (j){q, t^; cr) ~ fq, diverges at the glass transition as 
ta oc (T^* where 6 = l/2a. For short but macroscopic 
times, in the /3-regimc, one finds 

G^{t) cx {t/to)-", io < i < ta. (35a) 

For larger times, t > but still well below the so called 
a-relaxation time scale r 

G^(t) cx -(^/^^)^ U<t<^T. (35b) 
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FIG. 6. Complete set of exponents in 3D as function of the 
coefficient of restitution e. The top panel shows the exponents 
S and 7 for the divergence of the time scales. The lower panel 
shows the exponents a and b for the master functions and the 
exponent parameter A. 
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FIG. 7. Complete set of exponents in 2D as function of the 
coefficient of restitution e. The top panel shows the exponents 
S and 7 for the divergence of the time scales. The lower panel 
shows the exponents a and b for the master functions and the 
exponent parameter A. 



Finally, for the largest times the time-density superpo- 
sition principle holds, i.e., the coherent scattering func- 
tions can be collapsed on a master curve 

(f>{q, t; a) = 4>{q, t/T{q; a)). (36) 

The time scale r diverges at the glass transition, r = 
t((7;(t) oc (t^^ where 7 = (l/2a) -I- (1/26). An empirical 
choice for the scaling function (j){q, t) is provided by the 
Kohlrauseh- Williams- Watts law 0kww('J'jO ^ exp{—t^) 
[60] , crossing over to an exponential decay for the longest 
times [61]. 

The critical exponents a, b in Eqs. (35) are related to a 
single exponent parameter A = X{D,e) by the universal 
relation 

^ r^(i-a) ^ r^(i + 6) 

r(l-2a) r(l + 26)' ^ ' 

where T(x) is the Euler-Gamma- function. 

The sets of exponents shown in Figs. 6 and 7 in 3D 
and 2D, respectively, are functions of the dimension but 
differ only slightly between D = 2 and D = 3. Overall, 
there is a tendency for the exponents in Figs. 6 and 7 to 
show slightly lesser stretching for higher dissipation, i.e. 
smaller e. This may be interpreted as that the more dis- 
sipative and also more strongly driven fluid experiences 
less distinctive features in its glassy dynamics. The di- 
vergence in time scales is also expected to be a bit less 
pronounced. While the predicted changes in exponents 
are most likely hard to detect in experiments and simula- 
tion as absolute numbers, one should be able to detect the 
changes in comparison of different degrees of dissipation. 
Especially, the master functions should be comparably 
sensitive to changes in the exponents a and b. 



It is seen in Fig. 8 that for smaller e the order pa- 
rameter decays more slowly for large wave numbers 
indicating a tighter localization. In comparison to 3D, 
the transition in 2D exhibits individually sharper peaks 
and overall a tighter localization for the same dissipation, 
cf. Fig. 9. It has been shown for data from simulation 
[62] and experiments in colloidal suspensions [2, 63], that 
the MCT predictions for the are typically accurate 
to around 20%. Hence, rather than fitting individual 
directly to measurements and numerical calculations, 
experimental and simulation data can be expected to fol- 
low the distributions shown on the 20%-level and exhibit 
trends with variation of e as indicated here. 

It is seen in Fig. 4, that the critical law can only be 
observed without corrections for states closer than 0.1% 
to the transition point. Also, the von Schweidler law in 
Fig. 5 is only valid for an intermediate regime after the 
plateau. The regimes of applicability for the asymptotic 
scaling laws are therefore similar to the elastic case, and 
corrections to scaling are expected to follow the known 
trends [59]. 



V. TAGGED PARTICLE DYNAMICS 

The results of granular MCT for the tagged particle 
dynamics has been discussed in II. Here, we will focus on 
the derivation of the MCT equations. 

The incoherent scattering function, (pg {q, t) , captures 
the tagged particle dynamics. This includes the mean 
square displacement 5r'^{t) = (\jr s ~ r s{t)Y) which ap- 
pears as an expansion coefficient of the incoherent scat- 
tering function 0s(q, i) = 1 — g^(5r^(i)/6 -I- 0(q'*) and the 
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a macroscopic observable. 

We follow that reasoning and introduce the projector 



E 



(38) 



Together with the microscopic state a* = (Pq, jq^ /VT), 
it yields an equation of motion for 0s(<?j t), formally iden- 
tical to Eq. (16), 



ni 



dr ms{q, t 



Aq^T) 
^{q,T) 



(39) 



0, 



FIG. 8. Critical glass-form factors fg at the transition for 
e — 1.0 (dotted curve), 0.5 (full curve) and 0.0 (dashed curve) 
in 2D. 




with (j)siq,t = 0) = and 4's{q,t = 0) = 1. More- 
over, = {fq^\L+p'q) IVT = qVT and O^^- = qVT 
do not depend on the coefficient of restitution as shown 
in appendix B4. Hence, ft^ f2^„ilf,. = q^T is identi- 
cal to the corresponding quantity of the molecular fluid. 
This implies that looking at the probability density of the 
tagged particle on macroscopic time and length scales at 
very low densities, such that the memory kernel can be 
neglected, the microscopically broken time reversal sym- 
metry is unobservable. 

The memory kernels are given by 



iq,t) = (F^^(){t)F^)/q'T 



Uq,t)^{j^^U{t)F^)/qT. 



(40a) 
(40b) 



with the fluctuating forces = QgL+j^^ and F*^ = 
Qs^^jq^ and the fluctuating current is J^^ = Qs^^Pg. 



FIG. 9. Critical glass-form factors at the respective tran- 
sition points for e = 0.5 in 2D (full curve) and 3D (dashed 
curves) . 

diffusivity 6Doo ~ limt_>oo Sr'^{t)/t [51]. 



A. Equation of Motion 

Following the reasoning that one should account for 
the conserved quantities and only for the conserved quan- 
tities explicitly, one would assume that the equation of 
motion for a tagged particle should be flrst order in time. 
The density ps being the only conserved quantity as the 
momentum of the tagged particle is all but conserved. It 
has been shown, though, that a consistent treatment of 
the tagged particle dynamics in fact requires an equation 
of motion which is second order in time [64, 65], thus ef- 
fectively reintroducing the tagged particle momentum as 



B. The Mode Coupling Approximation 

We introduce a projection operator to describe the cou- 
pling between the tagged particle and the host fluid 

n^Nj2\pkPi){pkP-;\/Sk. m 

The corresponding mode coupling approximation reads 

k.p 

Within the MCA we flnd that L^iq, t) ^ and 
ms[(f>,(f)s]{q,t) w 

--^Y.^:ikpKkpmt)Mp,t), (43) 

k,p 
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where 

V^fcp = v/iV7^(jf |£+QsPfcP^) , (44a) 
W^fcp = v/iV7^(pfcP^|Qs£+J-f ) . (44b) 

Here, 

V^^^^ = ^={q-k){Sk-l)S^.k+p, (45a) 

is known from the htcrature [66, 67] [68] and one finds (cf. 
appendix B 6) 

Kkp = ^^^={q ■ k){Sk - l)Sg,k+p. (45b) 

For the vertices, the loss of detailed balance reappears 
also for the tagged particle dynamics. 
Together, Eqs. (45a) and (45b) yield 

m40, (j)s]{q,t) « 
1 + £ n f 

___ J d°kSk{q-kfcU{k,t)cj),{\q-klt). 

C. The approximate equation of motion 

Finally, the equation of motion for the incoherent scat- 
tering function reads 

'<Ps{q,t) +Vq(t)s{q,t) +nl(j>s{q,t) 

+ ni dTms[^,^s]{q,t-T)Mq,T)=0, (46) 
Jo 

capturing the coupling of the tagged particle dynamics 
to the dynamics of the host fluid as reported in II. 

The tagged particle is enslaved to the host fluid on 
macroscopic time scales. Consequently, it also develops 
persistent correlations, /| = limt_>oo 4>s{q,t), at the criti- 
cal density ipc- They can be calculated from the equation 

Y^^^lfJ-liq)- (47) 

Jq 

A detailed discussion of the solutions of Eqs. (46,47) was 
given in II. 

VI. DISCUSSION 

The granular MCT, which includes and extends MCT 
for elastic hard spheres, shows that the dynamics of a 
driven granular fluid is for one remarkably similar to the 
equilibrium dynamics and at the same time fundamen- 
tally different. It is similar in that there is always a glass 
transition, accompanied by the two step relaxation sce- 
nario of dynamic correlation functions and diverging time 
scales. As both the order parameter fq{e) and especially 



the critical exponents a(e) and 6(e) depend on the coef- 
flcient of restitution e, already slightly dissipative inter- 
actions (e ^ 1) destroy the universality of the dynamics 
on long time scales, which is observed in elastic systems 
with either Newtonian or Brownian dynamics [I, 69, 70]. 
The change of e will be detectable even on macroscopic 
time scales, in particular by observing the exponent b{e). 
This shows that the combination of dissipative collisions 
and driving cannot be mapped to an effective elastic hard 
sphere system with an effective temperature T^s, conceiv- 
ably different from the granular temperature, T. Such a 
mapping, would allow to find a scaling function such that 
(t){q,t]e) = 4>{q,t/toie)). 

The phase diagram in the (T, (y9)-plane is still an open 
problem. The jamming density is defined for athermal 
(T = 0) systems, while the granular glass transition is 
independent of temperature but assumes a finite temper- 
ature T > to sustain a fluid phase. It is not obvious, if 
and how they are connected. Our results suggest that the 
glass transition density, ipc, is always strictly smaller than 
the quasi static jamming density ipj. However, MCT is 
known to underestimate tpc- What happens for densities 
'P > fcis) larger than the critical density? At the glass 
transition, the a- relaxation rate diverges. Conse- 
quently, in every compression protocol using a small but 
flnite compression rate, at some density the compression 
rate will be larger than the a-relaxation rate t~^. From 
then on, the evolution of the system will be restricted to 
a subset of phase space. The packings, which are reached 
from that subset by further compression will also be re- 
stricted to a subset of all packings and that might not 
even include those of highest density. Even if the ideal 
glass transition is destroyed by processes which are ig- 
nored within MCT, the enormous increase of relaxation 
times will be prohibitive for all practical purpose. 

Apart from the Enskog term, Eqs. (19a, 19b), the equa- 
tions of motion are formally identical in two and three 
space dimensions. Hence the glass transition is qualita- 
tively similar in two and three dimensions, with however 
different values for the critical density and the critical 
exponents. Compared to = 3, the glass from factors, 
fq, decay slower in reciprocal space for D = 2, indicating 
a stronger localization in two space dimensions. 

In the flnal equations of motion, Eqs. (32) and (46), 
the driving force ^ appears only implicitly. While driv- 
ing is crucial to achieve a stationary state, beyond that 
it does not alter the relaxation rates flab, ^^'^ 'ioes 
not enter into the couplings to densities. Driving contri- 
butions would appear in the linear theory, if the (kinetic) 
energy, i.e., the granular temperature was included as a 
dynamic held. However this is hard to justify in a gran- 
ular fluid, where kinetic energy is dissipated locally and 
hence not a hydrodynamic variable. In terms of the MCT, 
a coupling to the currents in J'2 [Eq. (23)] would include 
explicit driving terms. Such a coupling was considered in 
the original mode coupling approaches [71, 72], but it's 
relevance even in equilibrium fluids remains unclear. 
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VII. CONCLUSION AND OUTLOOK 

We considered randomly driven inelastic smooth hard 
disks (in D = 2) and spheres (in D = 3). We systemati- 
cally derived equations of motion for the coherent scatter- 
ing function, (j){q,t), and the incoherent scattering func- 
tion, (l>s{q,t). 

The equations of motion are formally identical to the 
ones for elastic hard sphere or disk fluids in thermal equi- 
librium but acquire a nontrivial dependence on the co- 
efficient of restitution, e. A transition to a glassy state, 
indicated by a nonzero value of the order parameter, fq, 
appears through the bifurcation scenario of mode cou- 
pling theory. Like in thermal equilibrium, the spatial di- 
mension of the system only enters via the static structure 
factors. In both dimensions, the critical packing fraction 
increases the more dissipative the particles are. 

The dynamics around the plateau in the scattering 
functions is described by power laws with exponents, that 
are functions of the coefficient of restitution, e. Together 
with the e-dependence of the order parameter, fq, this 
shows that the dynamics fundamentally changes upon 
varying the coefficient of restitution. In contrast, the dif- 
ference between Newtonian and Brownian dynamics in 
thermal equilibrium, can be absorbed in the redefinition 
of the microscopic time scale. Also, a reduced long wave- 
length speed of sound is predicted for granular fluids. 

One can hardly expect to observe a glass transition in 
a fluid of monodisperse hard spheres, because the system 
would quickly crystallize. To slow down crystal nucle- 
ation, usually binary mixtures with a small size difference 
are used [73]. For fluids in thermal equilibrium, it was 
found that a MCT for mixtures does not yield results 
that differ drastically from those for the monodisperse 
idealization [74]. For mixtures of granular particles, a 
new complication will be the non-equipartition of energy 
between the mixture species [75, 76]. 

So far wc only derived equations for the correlation 
functions of spontaneous fluctuations. In a fluid in 
thermal equilibrium, this immediately entails knowledge 
about the corresponding response functions via the fluc- 
tuation dissipation theorem (FDT) [77]. In fact, a lot 
of the experimental measurements are concerned with 
response spectra [78] . The existence and form of a gener- 
alized FDT in driven granular fluids and more generally 
in systems far from equilibrium is a subject of active re- 
search [79]. 

It would certainly be desirable to weaken the assump- 
tions made on the stationary phase space distribution 
function g{T). So far, we ignore correlations between 
the velocities of different particles which are known to be 
present in driven granular fluids [80]. In light of the fact, 
that the single particle velocity distribution function is 
well represented by a simple Gaussian, these correlations 
can presumably be neglected as a first approximation. 
More serious are the static correlations, such as S{q), 
which are known from simulations to differ from their 
elastic counterparts used here. However, we can easily 



incorporate data for the simulated structure factors into 
our approach; work along these lines is in progress. 

The results presented above deal with a specific, highly 
idealized system. It is a natural question to ask how 
robust these results are qualitatively. No qualitative 
changes are expected for a speed dependent coefficient of 
restitution e = s{v). Also models that can be described 
by an effective coefficient of restitution [81] like, e.g., the 
spring-dashpot model are expected to show a nonequi- 
librium glass transition. The inclusion of inter-particle 
friction or the treatment of different driving forces will 
likely pose a number of challenges. Such changes might 
lead to equations of motion and results qualitatively dif- 
ferent from the ones discussed here. 
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Appendix A: The Granular Yvon-Born-Green 
Relation 



The Yvon-Born-Green (YBG) relation between the 
pair- and the triplet correlation function follows from the 
identity 



Vig(ri2) = {g ^Vi2£')5(7-i2) 



n j d^r3g3{ri,r2,r3){g ^Visp), 



(Al) 



where g ^ is the pseudo-inverse of the distribution func- 
tion [82]. For elastic hard spheres, where g(r) oc 



n.<, 



2a) This yields the known relation [49] 



Vi5(ri2) = ■ri2(5(ri2 - 2a)g{ri2) 

+ n rf^r3fi3(5(ri3 - 2a)g3(ri,r2,r3). 



(A2) 



For the inelastic hard spheres, there must be an addi- 
tional spatial dependence of the distribution function, de- 
pending on the coefficient of restitution, e, or otherwise, 
e.g., the structure factor Sq would not depend on e. 

Overlapping configurations still have zero probability 
and because of homogeneity, only relative distances play 
a role. Therefor, the distribution function will be of the 
form ^(r) cx rii<j 0(^u ^ 2a)'i?e(ry ) with a unknown 
hmction '!?£(r) > 0. With this, we get a granular hard 
sphere YBG relation 

Vig(ri2) = 1^12(5(^12 - 2a)5(ri2) -t- 5(^i2)Vi In-d^^r^) 
+ n (i^r3fi3(5(ri3 - 2a)g3(ri,r2,r3) 



<- J '^^'^3 53(^-1, r2,r3)Vi In i?e(ri3). 
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Unfortunately, virtually nothing is known about the use the approximation 
function i?e(r). Therefor we use the clastic hard sphere 
YBG relation which means we make the nontrivial ap- 
proximation 



g(ri2)Vilni?e(ri2) 

^-nj d^r3 53(ri,r2,r3)Vilni?,(ri3) (A3) 

which may be more general than setting '(9e(r) = 1. 
On the next level, wc have 

Vi53(123) = [fi2<5(ri2 - 2a) + f 136(^3 - 2a)]g3(123) 
+ g3(123)Vi[lnz?,(ri2) + ln,?,(ri3)] 



53(123)Vi[lnt9a^i2)+ln^?e(ri3)] 

~-nJ <i^r4 54(1234)Viln7?e(ri4) (A4) 

Appendix B: Matrix Elements 
1. The Frequency fijj in Two Dimensions 

We have to determine 



+ nj d^r4ri4(5(ri4 - 2a)g4(1234) 
+ n J d°ri 54(1234)Vilni?e(?-i4) 



where all other contributions vanish due to parity. Ex- 
using the abbreviation i = r^. In Eq. (B21) below, we plicitly, this reads 



J 



l + e 



(B2) 



where the three particle term vanishes, again, due to parity. Introducing the relative velocity v := {vi — V2)/\/2, the 
velocity averages can be evaluated 



{{q ■ f i)(fi2 • Vr2)'^Q{~ri2 ■ v^)) = V2{q ■ ^12) ((fi2 • f i2)^6(-ri2 • V12)) 



V2_ 
27rT 

-2rv/l>(g • fi2). 
I 



00 i'3it/2 

dv 



/2 



dip v'^ cos^ Lpe "^/^^ 



(B3) 



The remaining spatial average reads 



((q-ri2)^5(ri2-2a)(e*''-'---l)) 



((q-ri2)2<5(ri2-2a)(e^'J-'^- -1)) 

27r?iax 



N 



l + 24{2aq)]. (B6) 



2ax 



27r 



= cosV(l-e^""^'°"'^ 



) (B4) Collecting terms, one arrives at Eq. (19a). 

2. The Frequency Q,pj 



The driving contribution vanishes and the free stream- 
ing contribution yields 



One finds 



(B7) 



2ir 



dip cos tpe 



2 .^piz COS t/? 



. rfVo(^) 



= '?T5',/iV. 

(B5) The collisional contribution reads 
I 



N{N - 1) 



(Pq|7l2jq) 



1+f 



q ■ {{ri2 ■ ^^I2)'ri2e(-ri2 • ^;i2)<5(ri2 - 2a){e"^ ''' - e"^-^') ^ . e-'" ' 



(B8) 
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The velocity integration yields a factor T/2 while the 
spatial average can be rewritten as 



= 2(fi25(ri2-2a)e-*'J"--) 
+ 2{N - 2) (fi2(5(ri2 - 2a)e-''J"-^e"?-'-^) . (B9) 



Application of the YBG relation to the second term 
yields 



= ~^{S, - 1) - (fi2<5(ri2 - 2a)e-*«"--) , (BIO) 



i.e., the second term in Eq. (BIO) cancels the first term 
in Eq. (B9). Combining the remaining terms we arrive 



at Eq. (20b). 

3. The Vertex Wgkp 

As the vertex is linear in Vi, there is no contribution 
from the driving, i^C'jj. Expanding the projector Qc, the 
vertex reads 

Wgkp = N{pkPp\L+j^) /Sp 

- {pkpM {p,\^+3!^) /SpS,. (Bll) 

The free streaming contribution to the first term is given 
by 



(PkPpl-Cojq) = qT {pkPplPq} 



^5,,+p,qS^'Hk,p). 



(B12) 



For the coUisional contribution we find 



7V(^- 1) 



{PkPp\7t2J^) = ^^-^Tq ■ 



ik-Tj ^-ip rk 



ri2(5(ri2 - 2a) (e 



iq r2 piq r 







(B13) 



The average on the right hand side shall be abbreviated 
as (jfc|12). Then this can be expanded as 

(jfc|12) = (11|12) + (22|12) + (12|12) + (21|12) 
+ Ar((13|12) + (23|12) + (31|12) + (32|12)) 

+ 7V(33|12) + 7V2(34|12) (B14) 
Exploiting the symmetries, this can be simplified to 

(jA:|12) = 2 (11|12) + (12|12) + 2N (13|12) 

+ 7V(33|12) + 7V2(34|i2) (B15) 

Proceeding term by term we first find 

(11|12) = {Je'('=+P-'J'-'^ifi25(ri2 - 2a)(e*9"'i^ - l)^ 

which can be reduced to 

(11|12) = 4+p,q (ri2<5(ri2 - 2a)e*«"-i=) 

EE 5k+p,qG{q). 

The second term can be reduced to an equivalent cxpres- 



(B16) 



(12|12) = (e*^'=+P-'')"'iri2(5(ri2 -2a)e*'J-'' 

z'^''+P-''^-''^ri25{ri2 - 2a)e*('^-'»>''i^ 



(12|12) = 5k+p,q\G{k) + G{p)]. (B17) 
The first three particle term 



(13|12) ^e'('^+P-'')-''^e*'=-'^i^e*P-''^^fi2(5(ri2 -2a) 
_ (e'''k-i)-^^e'P-''Hi26{ri2 - 2a) J) 



r 



requires a little more work. The first term shall be ab- 
breviated as 



,i(;=+p-q)-r2gife-7-i2gip-r-32^^2(5(ri2 - 2a)^ 



H{k,p). (B18) 

The second term can be simplified with the help of the 
YBG relation 



e*('=-'j)''ie*P''^fi2(5(ri2-2a) 



—Sk+p,qMSp - 1) + NGip)]. (B19) 



Similarly, the second three particle term 
(33|12) = (|e'('=+P'-'^^ri2<5(ri2 - 2a)(e~'«"'^ - e~'i-ri^^ 
can be reduced by employing the YBG relation 



(33|12) - -—6k+pJiqiS, - 1) + NG{q)]. (B20) 

The four particle term 

(34|12) {e'''-'-'e'P-''^ri2S{ri2 ~ 2a){e-''' ''^ - e"'''"'!)) 

is naturally the most involved. Using the higher order 
YBG relation it reads 
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(34|12) 



2 
2 

iVT/3 



y d^r^d^'r-.d^u e-"J '-^e*'= '^^e^P '^V23^(r23 - 2a)53(r2, rg, r^) 
I d^r^d^r^d^'u e-*'J '-^e*^ '^^e^P '^V245(r24 - 2a)53(r2, rg, r4). 



Partial integration in the first term and the extraction of the momentum conservation constraint yields 
(34|12) = - ^4+p,q / d^r23d^r24 e-"= '-^^e-'P ''-.g3(r23, r-24) 



2 



Sk+p.g J d''r23d^r24 e-"= '-^^e-'P '--53(r23, r24)r23<5(r23 - 2a) 
Sk+p.g J d'^r^zd^'r^i e-''= '--e-'P '--.g3(r23, r24)r24'5(r24 - 2a). 



(B21) 



(B22) 



This leaves us with the relatively simple expression 



(34|12) 



2iq 
2 



N 



6k+p,q[H{k,p)+H{p,k)] 



l + e T 



[{q-k)Sp + {q-p)Sk-qS^^\k,p)]. 

(B24) 



Inserting Eq. (B12) and Eq. (B24) into Eq. (Bll) and 
applying the convolution approximation yields Eq. (29b) . 



Most terms cancel to yield 



2i 



4. The Frequency Q.i 



The free streaming contribution reads 



(pql^ojf ) = g((g-^s)')=9r (B25) 



(jA:|12) = —5k+p.q[kSp + pSk — qS^'^\k,p)] (B23) and the coUisional contribution 



{N - 1) (p^T+ ) = t^N{if,s ■ v,,fiq ■ risM-fis ■ v,s)6ins - 2a)) = 



(B26) 



vanishes due to symmetry. 



5. The Vertex V; 



qkp 



r 



as only the free streaming operator iLo applies. The 
velocity integration yield a factor of T 



The left incoherent vertex is given as 

Kfep = (jfl^+Pfc/^p) - (J'f l^+Pq) {P>kP;) ■ (B27) 
The triple density correlator, 



kT 



{jf\L+pkp-;) = — (p^Ip^P^) +pT(p^|pfcp;) 



N 
1 

N 



[{q ■ k)T + {q ■ p)Sk]Sk+p, 



(B30) 



(PqIPfcPp) = ■^Sk+p,qSk, 



(B28) 

is related to the structure factor. Moreover, we have 

(if |£+Pfcp;) = k {]f\:itPp) IPU'f ) (B29) Collecting terms one arrives at Eq. (45a). 
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6. The Vertex W^fep 
The incoherent vertex is given as 

W^fcp = {pkP'p\^+]t) - {PkPpK) (Pql^+jf ) • (B31) 
I 



can again be evaluated with the help of the YBG relation. 
Applying it to the second term cancels the first term and 

we get 

(iV-l)(pfcp^|T+jf ) 

= ^§4+p-q(g-fe)(5fe-l). (B35) 
Collecting terms one arrives at Eq. (45b). 

Appendix C: Details of the Numerics 

For the numerical solution of Eqs. (32,33,46), we used 
well established algorithms in 3D [59] and 2D [63]. Recip- 



The free streaming contribution is simple 

(Pfcp^l^ojf ) = qT (pfcP^jp^) = ^4+p,q5fe. (B32) 

For the coUisional part one finds with the velocity in- 
tegration being already performed 



(B33) 



(B34) 

I 

rocal space is discretized into grid points [M = 100) 
up to a cutoff of 2qa = 40 in 3D, and with M = 125 up 
a cutoff of 2qa ~ 50 in 2D. The time axis is also discrete 
with a grid of iV = 2048 points and a step size that is 
doubled in successive steps to accommodate for logarith- 
mic time scales. The initial time step is At = lO^^to. 
The critical density ipc is located by interval bisection. 



{N - 1) (PA,P^|T+ ) = (E, e-^^-^e-^^^-'^y^'iq ■ r.^ris - 2a)) 

The spatial average, 

(E, (9 • ^^Mris - 2a)) = (e-''^"-- {q ■ ns)S{ns - 2a)) + N (e'^'^-- {q ■ n^ris - 2a)) , 
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